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Abstract 

For most of the Ultra-High-Energy-Cosmic- Ray (UHECR) experiments and projects 
(HiRes, AUGER, TA, JEM-EUSO, TUS,...), the detection technique of Extensive 
Air Showers is based, at least, on the measurement of the air-fluorescence-induced 
signal. The knowledge of the Fluorescence-Light Yield (FLY) is of paramount im- 
portance for the UHECR energy reconstruction. The MACFLY experiment was 
designed to perform absolute measurements of the air FLY and to study its proper- 
ties. Here, we report the result of measurement of dry-air FLY induced by 50 GeV 
electromagnetic showers as a function of the shower age and as a function of the 
pressure. The experiment was performed at CERN using a SPS-electron-test-beam 
line. The result shows the air FLY is proportional to the energy deposited in air 
(Ed). The ratio FLY /Ed and its pressure dependence remain constant indepen- 
dently of shower age, and more generally, independently of the excitation source 
used (single-electron track or air shower). 
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1 Introduction 



The most challenging research in the field of cosmic-ray physics is certainly 
the highest-energy region of the cosmic-ray spectrum (around 10^'' eV) where 
the GZK effect ([l|) is predicted. The origin of these Ultra-High-Energy Cosmic 
Rays (UHECR) is still enigmatic. Precise measurements of their spectrum and 
arrival-direction anisotropy could provide important clue to understand their 
origin. However, the UHECR are rare, and their detection is an experimental 
challenge. In the 60's, a measurement technique based on the air-fluorescence 
light was proposed. This light, observed in the near-UV region (~300-400 nm), 
is induced by the de-excitation of the air molecules (mainly N2) occurring along 
the Extensive Air Showers (EAS) produced by the interaction of UHECR 
with the atmosphere. Since, most of the past, current and future experiments: 
HiRes (0), Pierre Auger Observatory (0), Telescope Array (4), Ashra (Q), 
OWL EUSO/JEM-EUSO (0), TUS Q), use the air-fluorescence signal 
to measure the UHECR flux, the knowledge of the Fluorescence-Light Yield 
(FLY) is of paramount importance for the UHECR energy reconstruction. 



In the past, the air fluorescence was a subject of extensive laboratory measure- 
ments. In 1967, A. N. Bunner summarized the existing data in his thesis (0) 
and proposed a FLY model with an uncertainty estimated at ~30 %. In spite 
of electron-beam-based measurements of Kakimoto et al. in 1996 (fl(5) and 



Nagano et al. in 2003 (jlll ). the uncertainties were still large inducing impor- 
tant systematics to UHECR experiments. A historical review of the air-FLY 
measurements can be found elsewhere dUl). 



The controversy and discrepancy between the AGASA (13 ) and HiRes (j2| 
experiments lead the community to pursue its effort to improve the knowledge 
of the air fluorescence (absolute light yield, spectrum) and of its dependencies 
with pressure, temperature, humidity, electron energy, shower age, etc. Since 
2002, a dozen of new independent experiments were proposed and took data. 
In these experiments, the air is excited with a high-energy-electron beam ( like 



for MACFLY, FLASH (jij) dlSl) and AirFLY (1161)), a radioactive source ([iTl) 
( flil ) (jigt ) or a low-energy-electron gun (20). A comparison and discussion of 
the results of these experiments can be found in the proceedings of the latest 
air-fluorescence workshop (21'). Two types of experiments were carried out, 
using thin or thick target. For the flrst type the fluorescence light is induced 
directly by interactions of primary electrons with the air, whereas for the 
second type, the light is induced by interactions of electromagnetic showers 



similar to EAS. Only two experiments used thick targets: FLASH (jl5l ) and 
MACFLY. 



The MACFLY (Measurement of Air Cherenkov and Fluorescence Light Yield) 
experiment has been designed to measure the light induced by both a single- 
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electron track and a high-energy electromagnetic shower developing in air (j22f ) . 
Actually, the experiment is composed of two devices, MFl for the single-track 
fluorescence (thin target) and MF2 for the electromagnetic-shower fluores- 
cence (thick target). The MFl results, obtained for different electron energies 
(1.5 MeV, 20 GeV and 50GeV), as well as an air-FLY model were previously 
published (1231). In this paper, we focus on the air-FLY measurements per- 
formed with the MF2 device at the CERN-SPS-X5-electron-test-beam line 
with 50 GeV electromagnetic showers. We compare these measurements with 
Monte Carlo simulations of shower development implementing a FLY model 
based on the MFl results. 



2 Experimental setup 



The MF2 device is composed of a pressurized chamber containing the gas 
under study, a pre-shower system and an optical system. The chamber is a 
quasi-cylindrical (96 cm in diameter and 146 cm long) large-volume (~lm^) 
tank internally covered with black paper (see Figure 1). The electron beam 
is aligned with the axial symmetry axis of the tank. It impinges on the pre- 
shower system which is a variable thickness target used to initiate electromag- 
netic showers inside the chamber. The pre-shower system is installed upstream 
the chamber, in a recess of 15 cm deep from the entrance wall of the tank 
(end-cap). The optical system measures the fluorescence light emitted by the 
excited air contained in the tank. It uses six UV-sensitive phototubes (PMTs) 
EMI9820QA. These PMTs are installed on the entrance end-cap on a 350 mm 
radius circle centered on the beam line. They point to the gas volume making 
a 20.5° angle with the beam direction. They are separated from the inner gas 
by quartz windows. The optical system is also composed of Winston cones 
and of different colored glass filters. The results presented in this paper cor- 
respond to the measurements performed using the Schott BG3 filters, which 
has a large transmittance band (290-440 nm). The same filter was used for the 
MFl measurements (The MFl device is described elsewhere ({23I) ) . 

The SPS-beam line is a pulsed beam delivering about 10 000 particles by spill 
(every 16.8 s). The purity of the beam is better than 98% for 50 GeV electrons. 
Along the beam line, upstream MF2, we have installed also the MFl device, a 
beam-position chamber (delay chamber) with a resolution of about 0.6 mm and 
a trigger system composed of two sets of scintillating counters, one upstream 
and one downstream the MFl device, as shown in Figure 2. The beam-spot 
size measured by the delay chamber is about 4 x 7mm^. 

An external gas system allows to fill both chambers, MFl and MF2, with the 
desired gas mixture, at a specific pressure. The two chambers are equipped 
with pressure and temperature gauges to control both parameters during the 
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PMT EMI 9820QA 




Fig. 1. The MF2 chamber: schematic view (left) and cut view (right). 




Fig. 2. Schematic view of the experimental devices along the test-beam line. 

data taking. In this paper, we report the FLY measurement of the following gas 
mixture: 80%(N2) and 20%(O2). This composition is close to the atmospheric 
dry air which is an admixture: 78.08%(N2)-0.93%(Ar)-20.99%(O2). The gas 
system enables to fill both chambers at the same time with the same gas, 
leading to simultaneous identical measurements with MFl and MF2. 



3 Pre-shower system 



When interacting in materials, high-energy electrons develop electromagnetic 
showers. The longitudinal development and the lateral spread of the shower 
are characterized respectively by the radiation length, Xq, and by the Moliere 
radius, Rmo, of the material (j24j ). In the air, an electromagnetic shower takes 
several kilometers for developing (at atmospheric pressure Xq — 300 m). Then, 
in order to sample air-shower development in laboratory, we have to use a fast 
and compact shower initiator. Here we use a pre-shower system made of a 
variable thickness copper target. 
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Fig. 3. The MF2 pre-shower system is a copper-disk-stack target surrounded by a 
lead-shielding tube. 

We choose copper because of its properties: high density (p = 8.96g/cm'^) and 
low atomic number (Z = 29) which initiate compact showers with a small lat- 
eral spread (Xo = 14.3mm, i?j\,/o = 14.9 mm). The characteristics of the elec- 
tromagnetic shower induced in the copper target are close to the EAS char- 
acteristics. The critical energy is of the same order of magnitude (~24MeV 
in copper and ~80MeV in air) and the energy of the secondary particles is 
similar. Moreover, the particle density of the shower in the MF2-chamber gas 
after the pre-shower (from 10^ to lO^e^/m^) is in the range of the particle 
density encountered in UHECR EAS at the shower maximum (from 10^ to 



Figure 3 shows a sketch of the pre-shower system. The target is made of a 
stack of copper disks, 10 mm thick each. The age of the shower changes as a 
function of the number N of disks in the stack. The equation 1 gives the pre- 
shower thickness (expressed in Xq units) as a function of N. One copper disk 
corresponds to (0.7 ± 0.002) Xq. All the matter on the beam line upstream 
the copper target (Trigger scintillator, MFl chamber, etc.) corresponds to 



In copper, at 50 GeV, the maximum of the shower development is at 7 Xq 
(~10cm). The pre-shower system allows to reproduce the air shower develop- 
ment in real atmosphere, on several kilometers, until the shower maximum. 
The sampling values used here are: 0, 1, 3, 5, 7 and 10 copper disks. 

In order to minimize the background induced in the PMTs by the showers, 
the copper- disk-stack target is surrounded by a lead shielding (20 mm thick) 



10^e±/m2) (j25t). 



(0.27 ± 0.05) X, 



Xn = (0.27 + iV X 0.7)X, 



(1) 
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which protects the PMTs from the backscattered particles. 



4 Data taking and FLY reconstruction 



The data recording is performed on an event by event basis. It uses a VME 
based DAQ system running a Labview program. The signal from all PMTs 
(MFl, MF2 and triggers) is recorded by QADC (CAEN-V792) which integrate 
the charge during 100 ns. We define two kinds of event: Beam Events (BE) and 
Random Event (RE). A BE is triggered when an electron passes through the 
two sets of scintillator counters (FLY measurement). The RE are randomly 
triggered (background measurement). For every run about one million events 
are recorded: 500 000 BE and 500 000 RE. 

The air FLY is rather weak and the majority of photons emitted are lost 
in the chamber. The typical mean number of photon detected by a PMT is 
about O.Olpe/evt (photoelectron per event). The method to extract the mean 
Detected Light (DL) of a run from the data is described in the MFl paper (f23|). 
This method can reconstruct the DL at the level as low as O.Olpe/evt with 
an uncertainty smaller than 4%. 

The detected light (DL) could come from several sources: Fluorescence (FDL), 
Cherenkov (CDL) and Background (Bgd). The overall signal is then: 

DL = FDL + CDL + Bgd . (2) 



Figure 4 shows the DL reconstructed from data and the estimation of CDL 
and Bgd contributions to the total measured light. The FDL is determined by 
subtracting CDL and Bgd to DL. The main part of the DL comes from the 
fluorescence whatever the conditions (50 GeV showers in dry air at 500 hPa 
for the left panel or at 100 hPa for the right one). 

The background level is determined from both, RE of the run and BE in 
vacuum where no light from neither fluorescence nor Cherenkov is expected. 
The background comes mainly from backscattered particles of the showers. 
The MF2 device has been designed to minimize it. As one can see in figure 4, 
the background is quite low compared to the fluorescence signal but it grows 
with the shower age (pre-shower thickness). That is why we limit our FLY 
measurements to the shower maximum. The uncertainty of the background 
measurement is about 20%. 



We estimate the Cherenkov radiation contribution with a Geant4 (1261 ) based 
Monte-Carlo simulation program. The Cherenkov light yield is important at 
atmospheric pressure (~ 20 ph/m/electron). However it is not contributing so 
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Fig. 4. Measured light in dry air at 500 hPa (left) and 100 hPa (right) in milli-pho- 
toelectron per event (mpe/evt) as a function of the pre-shower thickness (in Xq). 
Triangles represent the total signal (DL); dotted line is for the Bgd estimation from 
vacuum measurements; dot-dashed curve is the CDL simulation; stars are the FDL 
data (after substraction of Bgd and CDL); dashed line is the FLY model for showers. 
The solid line is the sum of the all contributions. 

much to the detected signal because it is mainly emitted in the forward direc- 
tion, downstream, where it is absorbed on the black surface of the chamber. 
The uncertainty of the CDL after diffusion on the black surface is estimated 
at 50%. 

Finally we extract the FLY in MF2 from the FDL, after dividing it by the 
MF2 efficiency eMF2 (see next section): 

FLY = D^-CI^^-Bod 

£mF2 



5 Calibration and systematic errors 



The Monte Carlo simulation plays a crucial role in the calibration and data 
analysis of the MACFLY experiment. The exact geometry and matter of all 
parts of MF2 (pre-shower, gas) and of other objects (MFl, scintillator) which 
are in the beam line have been reproduced carefully in a Geant4 based simu- 
lation program. For the external parts, only the details bigger than few cen- 
timeters are described. The optical properties of the optical system and of 
the inner chamber surfaces are also implemented. The fluorescence emission 
is assumed isotropic and its geometrical distribution is assumed to match the 
deposited energy distribution. The simulation program tracks all the optical 
photons until the PMT-photo-detection surface or until their absorption. 

All the phototubes used for the MACFLY experiment are tested and cross cali- 
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Errors sources 


Absolute 


relative 


MFl calibration 


16. 1 /o 




i\/rT7i /i\/rT7o 
iVll^ i / iVll^ Z 


io/o 




Geometrical distribution 


~ 0/0 


~ 0/0 


DL reconstruction 


~ 0/0 


~ 0/0 


CDL Simulation 


~ 1.5% 


~ 1.5% 


Bgd Measurement 


~ 1.5% 


~ 1.5% 


TOTAL 


23.1% 


~ 4.7% 



Table 1 

Systematic uncertainties of MF2 measurements in dry air at 100 hPa and for 5 Xq- 
thick-pre-shower target. 

brated in laboratory with a test bench using stabilized UV LED (370 nm) (22). 
Then, to calibrate the MF2 device we use the MFl chamber which is well cali- 



brated (1231 ). In order to do that, we performed measurement in both chambers 
filled with the same gas at same pressure and temperature, and without pre- 
shower target in front of MF2 (0 disk). In this configuration, both chambers 
measure FLY induced by the same 50 GeV-electron tracks. As the two mea- 
surements are in the same condition, the FLY/ should be the same in the 
two chambers. 



FLY/Ea 



FDLmfi 



FDLmf2 



£mfi X E, 



dMFl 



£MF2 X E, 



d.MF2 



(4) 



where FDLmfi and FDLmf2 are the fiuorescence light measured (in pe/evt) 
with MFl and MF2, Emfi and Smfi are the light collection efficiencies of 
MFl and MF2, E^^fi -E'dMF2 energy deposited in the air inside 

the measurement chambers of MFl and MF2. The collection efficiency of MFl 
is much better than the MF2 one. For the same light yield, the FDLmfi is 
about 12 times higher than FDLmf2- The ratio between the -EdA/Fi ^^"^ EdMF2 
is estimated by Monte Carlo simulation at 0.15. Then, we determine the MF2 
efficiency from the Emfi- 

eMF2 = ■ 1^ • eMFi = (7.0 ± 1.5)10-VMoton. (5) 

t JJLmFI ^dMF2 



The single-electron-track (no pre-shower) air fluorescence produces a very 
weak signal in the MF2 PMTs (< 0.001 pe/evt). The uncertainty on this mea- 
surement is large and induces an absolute calibration uncertainty of MF2 
larger than for MFl (see detail in table 1). 

The systematic error of -Edjv/Fi ^^'^ EdMF2i obtained by Monte Carlo simulation. 
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Dry air FLY: 50 GeV shower 
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Fig. 5. Dry-air-Fluorescence-Light Yield per unit of length (photon/meter) in dry 
air (100 hPa & 500 hPa) emitted by 50 GeV electromagnetic showers as a function 
of the shower age. Dotted lines correspond to a model of shower development in 
copper. 

is dominated by the air-density uncertainty (about 1% at 500 hPa and 2% at 
100 hPa). For the inter-calibration measurement, both chambers are filled with 
the same gas, then this systematic is negligible. However, the spacial geometry 
distribution of the energy deposited in the MF2 chamber changes as a function 
of the shower age. As the fluorescence emission should match the deposited 
energy distribution, the geometrical acceptance of MF2 could change as a 
function of the shower age. The systematic error induced by this effect has 
been estimated by simulation to be less than 3%. 

The background and the Cherenkov radiation represent a small fraction of 
the raw measured light (see figure 4) and induce small systematic errors. The 
Cherenkov contribution grows with pressure (density) and induces more un- 
certainty at higher pressure. Table 1 shows the contribution of the different 
systematic effects to the global uncertainty at 100 hPa and for 5Xo-thick pre- 
shower. The systematic errors of relative measurements are rather good (<5%) 
and grow only up to ~7.5% at 500 hPa. 



6 Result and discussion 

We measured FLY of dry air excited by electromagnetic showers for several 
pressures and for several shower-age values. Figure 5 shows the mean num- 
ber of fluorescence photons emitted when a 50 GeV shower traverses a one- 
meter-thick layer of air as a function of the shower age. We have performed 
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Fig. 6. Dry-air-Fluorescence-Light Yield per unit of deposited energy (photon per 
MeV) as a function of: (left) gas pressure, for two thicknesses of the pre-shower; 
(right) shower age, for two values of the pressure . A comparison with our FLY 
model (fi^ ) is also shown (dotted lines). 

such measurements for two different pressures: 100 hPa and 500 hPa. The dot- 
ted lines are proportional to a model of energy lost {dE/dX) by a 50 GeV- 
electron-induced shower developing in copper, based on Geant4 simulations. 
The measured air FLY follows well the expected shower development. 



To check the properties of air-fluorescence light induced by air shower, we 
compare our results to the air-fiuorescence model developed by Cohn (1221 ) 
found to reproduce well the MFl results (liil ). This model assumes the air 
FLY to be proportional to the energy deposited {Ed) in the air volume. Thus, 
the ratio FLY/ E^ (expressed in photons per MeV) should be independent of 
the excitation source. For each FLY measurements with MF2, E^ inside the 
chamber was obtained by Monte Carlo simulation. Figure 6 shows the varia- 
tions of FLY/ Ed as a function of the pressure (left panel) and as a function 
of the shower age (right panel). 



The pressure dependence was measured for two pre-shower thickness: 2.36 
and 5.16 Xq. We compare these results with the FLY model based on MFl 
measurements (1231 ) realized with the same gas mixture and the same optical 
filter (dotted line in Figure 6). For the two shower ages, the air FLY/ Ed 
has the same variation with pressure as for air fluorescence induced by a 
single-electron track. Comparison with other experiments is difficult because 
it requires a knowledge of the deposited energy in these experiments and 
corrections from the experimental differences (air composition, filter, etc.). 
However, the MACFLY model has been already compared with past (22) and 
recent (1211 ) experiments and shows similar pressure dependence. 



The shower-age dependence was measured at two pressures: 100 hPa and 
500 hPa. In both case, we do not find any significant variation of FLY/ E^ with 
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the shower age, in agreement with the FLY model (doted hnes). This result is 
also in good agreement with the FLASH thick-target experiments which mea- 



sured shower-age dependence in ambient air (jl5l ). This shows clearly that the 
FLY/ Ed properties are independent of the excitation source of the air. There 
is no clue of special behavior from saturation effect at the EAS density, or from 
low-energy-electron excitation of air molecules as it could be expected (j27|). 
Thus, air-FLY results obtained from electron-track experiments (thin target) 
can be directly used to determine the properties of the fluorescence induced 
by air showers. According to our result, the systematic error induced by the 
extrapolation from electron track to electromagnetic shower must be less than 
5%. 



7 Conclusion 



Using the MF2 device of the MACFLY experiment, we measured the Fluorescence- 
Light Yield induced by 50 GeV electromagnetic showers in a laboratory-controlled 
air (80%(N2)-20%(O2)). We studied both pressure and shower-age dependen- 
cies. The FLY variations with pressure, measured at two shower ages, are the 
same as the one measured with the single-electron-track device MFl (1231 ). The 
FLY variations with shower age, measured at two pressures, are well repro- 
duced by the shower-development simulations implementing our air-fiuorescence 
model which assumes the air FLY to be proportional to the energy deposited 
in the air volume. No evidence for FLY variation with the air-excitation source 
was found. 
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